ANALYTIC TORSIONS ON CONTACT MANIFOLDS 
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Abstract. We propose a definition for analytic torsion of the contact complex on contact 
manifolds. We show it coincides with Ray-Singer torsion on any 3-dimensional CR Seifert 
manifold equipped with a unitary representation. In this particular case we compute it and 
relate it to dynamical properties of the Reeb flow. In fact the whole spectral torsion function 
we consider may be interpreted on CR Seifert manifolds as a purely dynamical function 
through Selberg-type trace formulae. 



1. Introduction 



As introduced by Ray and Singer in [ 33 1 , the analytic torsion of a compact Riemannian man- 
ifold M may be seen as an infinite-dimensional analogue, on the de Rham complex (Q*M, d), 
of the Reidemeister-Franz torsion of a finite simplicial complex. More precisely, for A > 0, let 
E k be the ]0, A]-spectral space of Hodge-de Rham Laplacian A^ on /c-forms. Then the cut-off 
subcomplex (E* x ,d) is finite- dimensional and its Reidemeister-Franz torsion satisfies 

n n 

(1) 21nr R (Ei,d) = ln(jJdet(A klE ,)(-V k+lk ) = £(-l)*fcC'(A fcrf *)(0) , 

k=0 k=0 

where ((Ak^ E k)(s) = Tr(AfcjT^ fc ) is the truncated zeta function of A& on E\, Taking these 

(E^,d) as successive 'approximations' to the full de Rham complex, Ray-Singer defined the 
analytic torsion Trs as being 

1 n 

(2) r i?5 = exp(-^(-l) fc K / (A fe )(0)), 

k=0 

while the Ray-Singer metric on £ = det H*(£l*M, d) is given by 

(3) II \\rs = {TrsY 1 I k 2 (fi*M)> 

from the L 2 metric induced on C via identification of the cohomology by harmonic forms in 

n*M. 

The first purpose of this work is to adapt this idea to the contact complex (£* ,djf), a hypoel- 



liptic differential form complex naturally defined ([36, |37|]) on contact manifolds (M, H) of di- 
mension 2n+l. A specific feature of this complex is that the differential D = dn '■ £ n —> £ n+1 
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in 'middle degree' is a second-order operator, which is due to a slower spectral sequence con- 
vergence at this degree; see (3^, Proposition 3.3]. In order to find, as above, finite-dimensional 
cut-off subcomplexes (E k ,dn) approximating (£*,djj), we are led to consider fourth-order 
Laplacians in all degrees k; see (|lO|). The Reidemeister-Franz torsion of each cut-off 



subcomplex is then easily written (see Proposition [2.8D as 



2n+l 2n+l 



Ahxr R {El,d H ) = In J] det(A fcr£ ,) ( - 1) +1,u(fc) = £ (-l^MWH )' 



X 

k=0 k=0 



with w(k) = k for k < n and w{k) = k + 1 for k > n being the natural contact-weight of 
forms in £ k ; compare with (|l]). This leads us to define a candidate for the analytic torsion of 
the full contact complex by setting 

-, 2n+l 



k=0 



We define also a torsion function 

<s) = ^(-l) k+1 w(k)C(A k )(s) 



k=0 

and a Ray-Singer metric || \\c on the determinant of the cohomology det H*(£, dn), 

II lie = Tc | \l 2 {£) ■ 

(Our convention for Tc is inverse to Ray-Singer's original definition (§), but is standard now 
since it is natural at the metric level, as compared to (||); see e.g. |k| |TT||.) 

Having thus defined a torsion upon geometric and algebraic bases, we start then its analyti- 



cal study. We first establish in Theorem [3^ and Corollary general variational formulae for 
k(0), Tc = exp(K'(0)/2) and the Ray-Singer 'contact' metric || \\c- It turns out that k(0) 
is a contact invariant given by the integration of an unknown universal polynomial in local 
curvature data. Its vanishing is necessary in order for the Ray-Singer metric to be even scale 
invariant under change of contact form 9 i— ► K9 for K constant. We do not know whether 



k(0) vanishes in general except in dimension 3, as shown in Corollary |3.8| . Therefore the rest 
of the paper deals with this lowest-dimensional case. 



Corollary 3^ also states that there exist CR-invariant and contact-invariant 'corrections' to 
this Ray-Singer metric. Namely there exist universal constants (Cj)i<i<4 such that, on any 
contact manifold of dimension 3, 



| C R = exp(Ci/ R 2 0Ad9 + C 2 [ \A\ 2 9Ad9) 
\ Jm Jm J 



\c ■ 



is a CR-invariant (i.e. independent of contact form) metric on det H*(£, dn), where R and A 
are the Tanaka- Webster scalar curvature and torsion. Moreover 

|| ||£ = exp(CXM))|| ||cR and || ||g = exp(C 4 fftD*))|| IIcr 

are contact-invariant metrics, where v{M) is the ^-invariant of Biquard-Herzlich pf, an d 
rj(D*) is the CR-invariant correction to rj(D*)] see |7|, Theorem 9.4]. 
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We do not know the values of the constants Q. Though they are all related to the Heisenberg 
symbol of the Laplacians we consider, they might be difficult to compute: first since the 
are fourth-order, but also because the Heisenberg symbolic calculus (J5], 21], § |3.1|) suitable for 
the hypoelliptic contact complex, is highly non-commutative. 

Our next purpose is to compare the two analytic torsions and metrics coming from the de 
Rham and contact complexes. It is natural to expect they are related. Indeed, firstly these 
complexes have the same cohomology, being homotopy equivalent (Proposition |2.2j ), and in 
particular the determinants det H*(Tl*M, d) and det H*(£,dn) are canonically isomorphic. 
Moreover, from the point of view of spectral geometry, the non-exploding part of the Hodge- 
de Rham spectrum converges towards the contact complex spectrum when one takes the sub- 
Riemannian (or diabatic) limit e \ of calibrated metrics g £ = d0(-, J-) + e~ 1 9 2 ; see |38], 0|. 
Note that the classical Ray-Singer metric stays constant in this limit, being independent of 
the metric on M. 

However, we cannot prove equality of metrics in general, but only on particular contact 
manifolds called CR Seifert manifolds in (7j. These are CR manifolds (M,H, J) of dimension 
3 admitting a transverse locally free circle action preserving the CR structure (H,J); see 



Definition |4jJ . The generator T = d/dt of the circle action is the Reeb field of an invariant 



contact form 9. On such a manifold M, endowed with any unitary representation p : iri(M) 



U(N), Theorem L2 states that the two Ray-Singer metrics of the twisted de Rham and contact 
complexes coincide. 

In the last part of this work we analyse in detail the torsion function k(s) for CR Seifert 
manifolds. It first turns out that k(s) is a dynamical function in this case, depending only 
on the topology of M, together with the holonomies of the representation along the various 



primitive closed orbits of the circle action, as stated in formula (M) and Theorem 5.4. 



Specialising to s = leads in Theorem to an explicit formula for the Ray-Singer torsion 
and metric, twisted by any unitary representation. This Lefschetz-type formula extends a 
formula given by David Fried [[l?]] in the acyclic case, i.e. H*(M,p) = {0}, via topological 
methods and Reidemeister-Franz torsion. Fried interprets it as the identity 

(4) T RS (M,p) = |exp(Z F (0))|, 

where Zp(0) stands for the analytic continuation at s = of the dynamical function 



Z F (s) = -^ind(C Y )Tr(p(C))e- 



-sl(C) 

c 



Here the sum describes all free homotopical classes of closed orbits of the Reeb flow T, ind(C) 
denotes its Fuller index (Proposition |5.8|) , 1(C) its length and p(C) its holonomy. 



Our approach leads to another viewpoint on (|4j). Namely, we show in Theorem |5J| and (|66D 
that our purely spectral torsion function k(s) may be seen as a dynamical zeta function, in 
its whole. Indeed it holds that 

T( S ) (k(s) - k{M, p)) = - s)Z p {2s) , 
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for 

Z p ( S ) = ^ind(C) R Tr(p(C))£(C) s , 
c 

where again the sum runs over all free homotopical classes of closed orbits of the Reeb flow, 
R Tr is the real part of the trace, and n(M,p) = 2 dim H°(M, p) — dim H 1 (M, p) is a purely 
cohomological term. 

This Selberg-type trace formula also has a counterpart at the level of heat kernels. Indeed 
let Tr K (e~* A ) = 2Tr(e* A °) - Tr(e~* Al ); then we show in Theorem [5l0| that 



Tr K (e-* A ) = dimV^^ + -L V£(C) ind(C) R Tr(p(C))e~^ 2 / 4 * , 

where is the rational Euler class of the quotient surface orbifold E = M/S 1 . Hence our 
torsion heat trace (of fourth-order Laplacians) coincides with a dynamical theta function. This 
has some surprising consequences for the small time development of Tr K (e~* A ) on CR Seifert 



manifolds, but also on general 3-dimensional contact manifolds, as given in Corollaries 5.12 



and 5.13 



The paper is organised as follows. In §^| we first review one construction of the contact 
complex (£*,du) and recall the Ray-Singer argument, from the viewpoint of the Ray-Singer 
metric on the determinant of the cohomology. We then adapt this argument to the contact 
complex, which leads us to define the analytic torsion Tq, a torsion function k and a Ray- 
Singer metric || \\c on detH(£ *, du). 

In §|H we start the analytic study of this torsion. After reviewing relevant properties of 
hypoelliptic zeta functions and heat kernels, we establish variational formulae for re(0), the 
torsion Tq and the contact Ray-Singer metric || \\c- We then show that k(0) = in dimen- 
sion 3, and introduce corrections of the metric || \\c that give CR and contact invariants. 

In §[| we compare Ray-Singer analytic torsion to ours and show that the two Ray-Singer 
metrics coincide on CR Seifert manifolds. 

The final §|| is devoted to the study of the dynamical aspects of the torsion function of the 
contact complex, still on CR Seifert manifolds. 
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2. Contact analytic torsion via a determinant bundle 



Ray and Singer |33|| defined analytic torsion of the de Rham complex as an infinite- 
dimensional analogue of the Reidemeister-Franz torsion of finite simplicial complexes. Our 
purpose in this section is to describe their argument and adapt it to a similar complex defined 
on contact manifolds, the construction of which we now review. 

2.1. Contact complex. Let (M,H) be a smooth orientable contact manifold of dimension 
2n + 1. This means that the smooth contact distribution H C TM is given as the null 
space of a globally defined 1-form, called a contact form, satisfying the condition of maximal 
non-integrability 9 A {d0) n ^ 0. The contact forms comprise an equivalence class under 
multiplication by smooth non-vanishing functions. 

The contact complex (|36], ^]) is a refinement of the de Rham complex on contact manifolds 
defined as follows. Let Q*M denote sections of the graded bundle of smooth differential forms 
on M, I the ideal in Q*M generated by 9 and dO, and J the ideal in f2*M consisting of 
elements annihilated by 9 and d9. One verifies that Z k = £l k M for k > n + 1, J k = for 
k < n, and that the de Rham exterior derivative d naturally induces operators dn to form 
two complexes 

n°M ^ nHi/i 1 Q n M/l n 

and 

j-n+l d H ) j-n+2 _ _ _ d H ) j-2n+l 

It is clear that these two complexes are defined independently of the choice of 9. These two 
complexes are joined by a second-order differential operator D : Q n M/T n — > J n+1 defined by 
setting D[a] = d(3, where (3 £ S1"M is defined by the following: 



Lemma 2.1 (|3(| |37j). Let a G Vt n M. Then there exists a unique (3 G £l n M such that (3 = a 



(mod 9) and 9 Ad/3 = 0. Moreover d(3 G J n+1 , and if a G l n then dp = 0. 

One can show that D may in fact be defined independently of the choice of 9. The contact 
complex is 

n°M ^ rtM/i 1 ^ ... ^ n n u/i n j n+1 ^ j n+2 ... ^ jzn+\ 

We will write £ k to mean either £l k M/2 k or J k . We also have: 

Proposition 2.2 (|37], pg. 286]). The contact complex forms a resolution of the constant 
sheaf M. and hence its cohomology coincides with the de Rham cohomology of M . Moreover 
the canonical projection ir : Q k M — > Q k M/T k for k < n and injection i : J k — > ft k M for 
k > n + 1 induce an isomorphism between the two cohomologies. 

The arguments being purely local, these results also apply to twisted versions of the complex 
with a flat bundle, as coming from a representation p : tti(M) — > U(N). 

It is a basic fact that the symplectic bundle (H, d9) admits a contractible homotopy class 
of calibrated almost complex structures, i.e. J G End(H) is in this class if and only if J 2 = —1 
and the Levi metric d9(-, J-) is positive definite and Hermitian. 
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The Reeb field of 6 is the unique vector field T satisfying 9{T) = 1 and T_i d9 = 0. Fixing 
a 6 and a J, we may define a Riemannian metric g on M by using the Levi metric on vectors 
in H and declaring that the Reeb field T is of unit-length and orthogonal to H, i.e. 

g = d8{-, J-) + 6 2 . 

With these choices, it is straightforward to identify the quotients of forms appearing in the 



lower-half of the contact complex with actual forms (||37|, pg. 288]): 

n k N/l k = {a G n k N : Tj a = = Aa}, 

where A is the adjoint of the operator L : Q k N — > fl k+2 N, La = d6 A a. 

We henceforth assume that M is compact. With the identification above we now have 
an L 2 inner product defined on the contact complex. Let djj, D* denote the formal adjoint 
operators. It is straightforward to verify that 

(5) 8 H \ £ k = {-if *d H *, D* = (— l) n+1 * D*, 

where * : £ k —> £ 2n + 1 ~ k [ s the usual Hodge * operator. 

As a last comment here, we mention there exist other approaches to this elementary con- 
struction of the contact complex. One possibility is via spectral sequence considerations, using 
a canonical filtration by Heisenberg weight of forms 0*M; see |38], §3] and 0| where this ap- 
proach is used in the study of the sub-Riemannian (diabatic) limit of the Hodge-de Rham 
spectrum. Another interesting viewpoint is to consider the contact complex as a curved ver- 
sion of a Bernstein-Gelfand-Gelfand complex in parabolic geometry; see |2[ §8.1] for such a 
presentation on the 3-dimensional Heisenberg group. 

2.2. Determinant bundles, metrics and Reidemeister-Franz torsion. We follow the 



presentation of Bismut and Zhang 1 13] to define the Reidemeister-Franz torsion of a finite- 
dimensional complex. 

Let £ be a finite-dimensional real vector space, and define the line 

det£ = A max £. 

A useful convention here is to set det{0} = R (compatible with det(i? © F) = det E (g> det F). 
If A is a line, let A -1 = A* be its dual line. Then A (g> A -1 = End(A) = Rid is canonically 
isomorphic to R. 

One extends these notions to a finite-dimensional complex. Let 

(E,d) : —> E Ex • • • E n — 
be such a complex and H*(E,d) its cohomology. Define 

n 

detE = (g)(det£ fc ) ( - 1)?c 

fc=0 

and 

n 
fc=0 



det{H*(E,d)) = (g) (det H k (E,d))( 
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Proposition 2.3 (Knudsen-Mumford |2q1 ). The lines detE and det(H* (E , d)) are canoni- 
cally isomorphic. 

Proof. We include the proof as the explicit form of the isomorphism will be useful below. We 
follow jT(|- Suppose first that H*(E,d) = {0} so that det H*(E, d) = R. Then we need to 
find a canonical section of detE. Let = dim Ef.. Pick a non- vanishing element 

(6) s = ei A e 2 A • • • A e No G A^^o = det E ; 
then 

(7) ds = dei Ade 2 A ••• Ade No G A^-Ei 

is non- vanishing since d : Eq — > i?i is injective. 

Next pick si G f\ Nl ^ N °E\ such that dso A s\ generates deti?i, and so on, taking Sk G 
^iv 4 f-(-i) N o]^ i such that dsk-i A generates deti^. Consider now 

(8) S(E, d) = s <g) (ds A si) -1 <g) (dsi A s 2 ) (8) • • • <8> (ds n _i) ( ~ 1)n G det E. 

It is clear that the class S(E,d) is non-zero and does not depend on the choices of Sk for 
k = 0, . . . , n, completing the proof of the proposition in the acyclic case. 

For the general case, observe that the determinants of the short exact sequences 

-»■ dE k -» kerd| £ ; (s+1 -> H k+1 (E,d) -► 
-> kerd|£; fc+1 -> £fc + i -> ->• 

each have a canonical element, as was just shown above. So we have canonical isomorphisms 

det(kerd|£ fc+1 ) det(dE k ) ® det(# fc+1 (£, d)) 
det(£ fc+ i) det(kerd| i3fe+1 ) ® det(dE fc+1 ), 

and then 

det( J B A;+ i) det(dE k ) ® det(H k+1 (E, d)) <g> det(dE fe+ i). 
Finally taking tensor products over k gives 

det So ® (det iii)" 1 <g) det E 2 ® • • • = det 77° ® (deti? 1 )" 1 (g> det H 2 ® ■ ■ ■ 
canonically. □ 

Suppose now i? is given a metric g. Hence det E has an induced metric. One can then 
define a metric on det H*(E,d) by 

I Wdet H*(E,d) = II l|det£> 



using the canonical isomorphism given by Proposition 



Let d* = 5 be the adjoint of d. By finite- dimensional Hodge theory, H*(E, d) identifies with 
the harmonic forms 

H*(E,d) = {se E\ds = d*s = 0}. 

By their inclusion in E, the harmonic forms inherit a metric. We then have a a second metric 
\detH*(E,d) on det H*(E,d) via the above identification. 
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Definition 2.4. The torsion of the complex (E,d) with metric g is the ratio 

z -w—! j \ I! Ildet H*(E,d) 

t(E, d, g) = - — . '- . 

I Idet H*(E,d) 

Remark 2.5. Note that this definition of torsion, given in [l3|, §2] for instance, is quite natural 
at the metric level, but actually leads to the inverse of the original Reidemeister-Pranz torsion 
(or .R-torsion), i.e. 

r R (E,d,g) = l/r(E,d,g); 

see [||, §1] and 0, §2]. 



rMufMr 1 x---xf„ IK - ill ,.V' ir 

as n -i\ 



One can be more explicit using the proof of Proposition |2.3| . Consider 

F = H*(E,d) ± . 

The complex (F, d) is acyclic so we can construct the canonical class S(F, d) as in 
Proposition 2.6. For any generators S{ E det(kerd) -1 " we have 

(9 ) TW4»>HsW«0l-*-(g)x(g) 

Proof. The splitting E = H.*(E,d) © F induces the canonical isomorphism 

det H* (E, d) det £ = det H* (E,d)® det F 
si — ► s ® S(F, d) . 

Then 

\\s (g> S(F,d)\\ detE = \s\detH*(E,d)\\S{F, d) Ildet F 

and by Definition 

T(E,d,g) = \\S(F,d)\\ dctF 

= \\sq\\ x \\dso A si|| _1 x \\dsi A Ss || X • • • X ||( _1 ) 



n-l 



= (JM.) x f JfiLr 1 x ... x f l|Sn ~ l11 V' 1] 

V ||c?so II ^ Ml^ s ilr Ml^ s n-ilr 
since if S{ S det(ker d)- 1 , \\dsi A Sj+i|| = ||(isj||||sj + i||. □ 



At this point in the Riemannian case ( |p3[ |13|| ) one can guess the correct formula for analytic 
torsion by considering the Reidemeister-Franz torsion of finite-dimensional subcomplexes that 
approximate the infinite-dimensional de Rham complex (Q*,d). A natural choice of subcom- 
plexes is obtained here by taking cut-off de Rham complexes using the spectrum of the Hodge- 
de Rham Laplacian A = dd + 5d; that is, one considers the energy levels A] = — 
One then expresses (§) using the determinants of A on {SV^^d) and finally as combinations 
of differentiated zeta functions (' (A)(0) in the limiting infinite-dimensional case. We carry 
out this procedure for the contact complex next. 
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2.3. Defining a contact analytic torsion. Consider now the contact complex 

(£ djj) • £ dH > E 1 da > ■ ■ ■ dn > £ n D > £ n+1 d{1 > ■ ■ ■ dn > g 2n+1 

We want to define finite-dimensional subcomplexes of the contact complex via finite energy 
cut-offs for a certain Laplacian A. 

We shall use the following uniformly fourth-order Laplacian: 



(10) A ■, 



{d-H^H + 5h(Ih) 2 on £ k for k / n, n + 1 
{d H 5 H ) 2 + D*D on£ n 
DD* + (5 H d H ) 2 on£ n+1 . 



We denote by A^ the restriction of A to £ k . The rationale behind our choice of A is as 
follows. In middle degrees, because D is second-order, one needs to square the terms involving 
du so that A has certain good analytical properties. In particular, A is maximally hypoelliptic 
and invertible in the Heisenberg symbolic calculus, while the standard combinations d}jo~H + 
D*D and DD* + dn^H are not; see § |3.1| below. If we then consider the spectral spaces 
EL, = {A n < A} as successive finite-dimensional approximations of £ n , in order to include 
E™ Q A j in a finite-dimensional subcomplex of (£*,d) we need the Laplacians outside middle 
degree to be fourth-order also. 



Remark 2.7. Note that A is different to the Laplacian Aq defined in p7|| . The latter was 
defined with nice algebraic properties, namely commutativity with J when Ct J = 0. On the 
other hand, observe that A commutes with dn,D and their adjoints, which, as we shall see 
in § [3^ , is essential for analytic torsion of the contact complex having the correct variational 
behaviour. Moreover note that A is the Laplacian appearing in the sub-Riemannian limit 



Next let us set 

2n+l 

E [0M= 0{A fc <A}. 

fc=0 

These are finite- dimensional subcomplexes of the contact complex, each of which splits into a 
direct sum of eigenspaces E*^ of A for < fi < A. 

Proposition 2.8. The torsion of (Ef Q x\^h) is 

2n+l 

(11) r(E^ X] ,d H ) = J] det(A fe | E^ x] )^ k ^/\ 

k=0 

where 

k if k < n 



(12) w(k) 



k + 1 if k > n . 



Remark 2.9. This w(k) is the natural weight of £ k in the contact complex; see §3J_ and 
§3]. 
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Proof. We use Proposition 2i5. By the orthogonal splitting EL^i = (B^<\E^ and @ we have 

(13) T (E^ x] ,d H )= n 11^,^)11. 

0<m<A 

Given [i G]0, A] we pick an orthogonal basis (t>i, • • • , U7v fe ) of each eigenspace 

^n(kerd) ± = F / ^ 

and choose the elements s\ = v i A • • • A VN k G det as in d). Here iVfc = dimF*. Now by 
(0) and orthogonality 

\\dsi\\ = \\dv\ A dv2 A • • • di>7v fc || 

= \\dvi\\ x 1 1 rfi?2 1 1 x • • • x ||dvjv fe || 

{^A r fc/4|| Ul || x ... x \\vN k || = fi Nk 7/4 1 1 Si|| for k ^ n 
[i Nk l 2 1 1 ui|| x • • • x ||t!Ar fc || = [i Nk l 2 1 1 Sj 1 1 for k = n . 

Therefore from (g) 

(14) \\S{E^d H )\\=^ 
with 

4M„ = -N + iVi - iV 2 • • • + (-l) n+1 2iV n + (-l) n iV n+ i ■•■-N 2n . 
Let ek = dim£^. From = © dF^ 1 , one has = JVjfe + Nk-i and therefore 



This leads to 



N k = e k - e fc _i + ■■■ + (-1)%. 

4M M = -(2n + 2)e + (2n + l)ej + • • • + (-l) n (n + 3)e n _i 
+ (-l)" +1 (n + 2)e n + (-l) n+2 ne n+1 
+ (-l) n+3 (n - l)e n+2 + • • • + e 2n+ i 

2n+l 

= ^(-l) fc u;(A;) efc , 

fc=0 



2n+l 



using that (— l) fc efc = for the acyclic complex (E*, da)- 



As A = n on E* equation (O) reads 



2n+l 



||S(i^)|| 4 = H det(A fc I ^) (-1) 



fc=0 



leading by ( jl3|) to the assertion of the proposition. □ 

We finally introduce zeta functions of the contact Laplacian. If spec*(A*.) denotes the 
non-zero spectrum of A& on then we take 

C(A fe )(s) = dimtf fc (£, d H ) + X ~ S - 

Aespec*(A fe ) 
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Note that by hypoellipticity (or Proposition |2.2| ) dim H k (£ , dff) is finite. By the results in 
§ 3.1 below, ((Ak)(s) admits a meromorphic extension to C that is regular at s = 0. On each 
subcomplex (E* Q ^,dn) we then have 

C'(A fe |£ ]0 , A] )(0) = - Yl ln v 

^espec*(A fc )n]0,A] 

= -lndet(A fe | £] ,a]) . 
Thus formula ( [TT| ) for the torsion of (E* Q X ]^h) can be written for A > as 

-, 2n+l 

(15) lnr(E^ x] ,d H ) = - £ (-l^V*)^* | £f 0>A] )(0) . 

fc=0 

We thus speculate in extending this formula to the whole contact complex by defining the 
analytic torsion of the contact complex as 

, 2n+l 

(16) lnT c = I ^(-l) fe+1 ^)C'(A fc )(0). 

fc=0 

This formula is very similar to that of Ray-Singer analytic torsion Trs in the Riemannian 
setting. Namely, from [ 33 , Definition 1.6], in dimension N 

1 - 

(17) lnr RS = -^(-l)^C'(A fc )(0), 

k=0 

for Hodge-de Rham Laplacians A&. Note however the sign convention: Trs coincides with 
Reidemeister-Franz torsion tr on finite- dimensional cut-off de Rham complexes, while our Tc 
leads to the inverse; see Remark 2.5. 



By analogy with Definition 12^ and |32], 13|, we also define a contact complex Ray-Singer 
metric on detH*{£,dn) by setting 

(18) II lie = Tc | \l 2 {£) ■ 

Here | \l 2 {s) 1S the L 2 metric induced on det H*(£, dn) by identification of H*(£,dn) with 
harmonic forms H.*(£,di{) C £* . 

Again, note that the Ray-Singer metric on the de Rham determinant det H*(Q*M, d) reads 
instead 

(19) II Wrs = {Trs)~ 1 \ Il 2 (q*a/) • 

More generally, we can twist the contact complex with a flat bundle and then define the 
analytic torsion of this twisted contact complex (£*,dn). Indeed let p : iri(M) — » U(N) be a 
unitary representation on C N . Associated to p is an Hermitian complex rank iV vector bundle 
V p equipped with a canonical metric-preserving flat connection V p . One sets £ p = £ ®V p and 
dn(a s) = dna (g) s for parallel s. From this we may define the contact analytic torsion 
Tc{p) with associated contact complex Ray-Singer metric on &etH*{£p,d}{)- 

The conciseness of notations Tc and Tc{p) should not be misleading. The (twisted) contact 
complex only depends on the contact structure H on M (and p), but the spectral invariants 
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Tq and Tq(p) also depend on the choices of a contact form 9 and complex structure J, both 
being used in the metric g. 

Since we have defined this analytic torsion through algebraic and formal considerations 
around Reidemeister-Franz torsion, we now need to study its analytical properties. That is 
the purpose of the next section. 



3. Heat kernels and variational behaviour of the torsion 

We first gather some properties of zeta functions and the heat development of hypoelliptic 
operators such as the Laplacian of the contact complex. 

3.1. Heat kernels and zeta functions for hypoelliptic operators. The Laplacian A for 
the contact complex is not elliptic. However there is a (substantially more intricate) symbolic 
calculus that can be applied to it to obtain results on heat kernels qualitatively analogous 
to the elliptic case. This calculus is called the Heisenberg calculus and was introduced by 



Beals-Greiner [|| and Taylor [44]. A short account of its properties may be found in pl}| , and 
its use for the contact complex has been presented by Julg and Kasparov in |^3|, §5]. This 
calculus has also been developed in a more general setting that includes the contact case; see 
the presentation of Ponge (3(| for details and further references to the literature. Here we just 
briefly sketch the results that we shall need in the sequel. 



Proposition 3.1 (Ponge |30|, Ch. 5]). LetV be a vector bundle over a compact contact manifold 
(M, H) of dimension 2n + 1. Let P : C°°(M, V) -> C°°(M, V) be a differential operator of 
even Heisenberg order v that is self adjoint and bounded from below. If P satisfies the Rockland 
condition at every point then the principal symbol of P+dt is an invertible Volterra-Heisenberg 
symbol and as t \ the heat kernel kt(x, x) of P on the diagonal has the following asymptotics 
in C°°(M, (End V) ® |A|(M)): 

oo 

E 2(j — n — 1) 
t v dj(P)(x). 

3=0 

Some explanation about the proposition is in order. The Heisenberg order of P is defined by 
assuming that a derivative in the direction of the Reeb field T has weight 2, while derivatives 
in the direction of the contact distribution H have weight 1. The Rockland condition is a 
representation-theoretic condition defined in Definition 3.3.8]. (The original formulation 
is due to Rockland p4|.) An operator that satisfies this condition is hypoelliptic, in the sense 



of [30, Proposition 3.3.2]. Invertibility of an operator in the Volterra-Heisenberg calculus is 



explained in [30, Ch. 5]. 



The next result describes the properties of the zeta function in this contact setting. For the 



non-negative operators P we are concerned with, the result follows from Proposition |3.1| by a 
classical argument using the Mellin transform of the heat kernel. 



Proposition 3.2 (Ponge |3T|, §4]). Let P be as in Proposition \3.l\ . Then the zeta function 



C(P)(s) = dimkerP + Tr*(p- s ) , s G C, 



ANALYTIC TORSIONS ON CONTACT MANIFOLDS 



13 



is a well-defined holomorphic function for Re(s) S> 1 and admits a meromorphic extension to 
C with at worst simple poles occurring atseS = { 2(n+ t) 1 ~ j) | j e N}\ (-N). Moreover 

C(P)(0)= / tr( a „ + i(P))0A(d0) n 
JM 

is the constant term in the development of Tr(e~ tp ) as t \ 0. 



Now by |37|, p. 300] and the definition (|30|, Definition 3.3.8]) of the Rockland condition, 
or else by Julg and Kasparov's work |23|, §5], the fourth-order Laplacian A on the contact 



complex (twisted with a flat bundle) satisfies the Rockland condition, hence Propositions 3J. 
and |3^ apply to it. 



Remark 3.3. Our convention is to include the eigenvalue in the definition of the zeta function 
of a non-negative hypoelliptic operator P. This means C(P)(0) consists of (the integral of) 
purely local terms for P. 

3.2. Variational behaviour of the analytic torsion. We consider the variation of contact 
analytic torsion and Ray-Singer metric in the direction of an arbitrary line of pairs (9 e , J £ ) 
of contact form and calibrated almost complex structure for (M,H). We shall see that the 
variation of the Ray-Singer metric is given entirely by local terms. Indeed this may be 
viewed as a necessary and sufficient condition for correctly defining an analytic torsion; see 
the approach of Branson [l4|] . 

First we define from ([if]) the contact torsion function by 

2n+l 

(20) «(*) = - ^(-l) fe+1 ^)C(A fc )(,), 

k=0 

with w(k) as in (12). Then the analytic torsion of the contact complex reads 

lb = exp(-«'(0) 

For simplicity, in this section we suppress the C from the notation, as well as the represen- 
tation p, although all results stand for the twisted torsions and metrics as well. 

Theorem 3.4. Let a • superscript denote first variation (d/de)\ £= o. 

(1) One has n(0)' = 0, so that k(0) is a contact invariant. 

(2) The variation of the analytic torsion T is given by 

(InT)' = V(-l) fc ( [ tr(aa t0 , k ) dvol - Tr(aV k )) , 

where a = a t o. k is the t° coefficient in the diagonal small-time asymptotic 

expansion of the heat kernel of A/%, dvol is the volume form A {d0) n , and V k is 
orthogonal projection onto the null-space of A k . 

Proof. By Hodge * duality (for convenience we suppress the e dependence) 

n 

K ( s ) = ^2c k C{A k )(s) 

k=0 
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c k = (-l) k (n+l-k). 



By a Mellin transform set 



71 /*+oo n 

f(s) = T(s)K(s) = Y J Ck / t^Tiie-^ -V k )dt + T(s)Y,Ck dim H k (£,d H ), 

k=0 ^° fc=0 

where V k denotes orthogonal projection onto the null-space of A&. 
Lemma 3.5. It holds that 

(Tr(e-' A ))' = -tTV(A*e-* A ). 
Proof. Duhamel's formula (see e.g. |3f| Proposition 3.15]) implies that 

Tr(e-*( A+£A *+°( £2 ))) - Tr(e"* A ) = - /"* Tr ( e -(*-s)(A+ e A-+o( £ 2 )) e(A . + (e)) e - sA ) ds. 

Jo 

Dividing both sides by e then letting e — > gives 

(Tr(e-* A ))* = - I' TT(e-^- s)A A'e- sA ) ds. 
Jo 

Recalling that Tr(AB) = Tr(BA), for smoothing operators A,B, we have 

(Tr(e- tA ))* = - /** Tir(A'e- aA e-(*- a ) A ) ds, 
Jo 

and the lemma follows. 



□ 



By Lemma p.5 

n r+oo 

(2i) /(*)' = -£<*/ * STV ( A ' f 

fc=o Jo 



-tA k 



since by hypoellipticity Tr(V k ) = dimker A^ = dimiT (£, du) is certainly independent of 6 
and J. 

Setting a = one computes using @ the variation of the Laplacian as 

—dHCt^HdH^H + dndHOedHdH — dH^HduO-^H + du^HdH^H^ 

—adHdnSndH + ^HOtdn^HdH — SjjdHCtSjjdjj + dudu^H^dH 
on £ fc for /c = 0, . . . , n — 1 , 



A* 



-dHOidudH^H + du5}jOid}j5}j — dHSndHO.dH + dH^Hdn^H^ 
—aD*D + D*aD 
on £ n . 



A computation then shows that 



ra-1 



(22) 



]T Cfc Tr(A£e-* Afc ) = 2^Tr {{a(c k + c k ^)(d5) 2 - a(c k + c fc+1 )(<M) 2 ) e ~ tA *) 



fc=0 



fc=0 



+ 2Tr ((a(c n + c n _i)(^) 2 - ac n D*D) e~ iA ") 
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To move the a's to the front we have used the following facts: the heat kernel is a semigroup, 
implying e~ tA = e ~ (t/2)A e~ ( * /2)A ; if operators A,B are smoothing then Tr(AB) = Tr(BA); 
and the Laplacian A commutes with dn,D and their adjoints. We have also used that 



Ti(aDD*e- tAn+1 ) = -Tr(<xD*De 



as aDD*e * A ™+ 1 = - * 1 (aD*De * A ™)*, which follows from (|) and (* 2 )* = *a + a* = 0. 
Simplifying (53) yields 



^c fc Tr(A^e-' Afe ) =2^(-l) fc+1 Tr(aA fc e- 



k=0 fc=0 



= - 2 |E(- 1 ) fe+lTr (« e " tAfc )- 

fc=0 

Hence after integrating by parts in @), 

n /•+oo 

(23) /(a)' = 2s V(-l) fe / t s_1 Tr(ae"' Afc - aV k )dt . 

We have included the aVk term in the integrand to ensure that f(s)' is a well-defined mero- 
morphic function on all C. 

As l/r(s) = O(s), f(s)' is regular at the origin. Now near s = 

/( s ) = r( s )( K (o) + SK '(o) + o( s 2 )) 

so 

f( s y = r(s) K (py + s r( s )(«'(o))* + o( s ). 

Since f(s)' is regular at the origin, letting s — > in the above equation gives 

k(0)* = and (k'(0))* = /(0)V 



This gives assertion (1) of Theorem |3.4| , Assertion (2) follows from Proposition |3.2| , after 



taking the inverse Mellin transform of ( p3| ) . □ 
Remark 3.6. The previous proof actually shows that the torsion function 



K 

k=0 

studied here is, up to a multiplicative factor, the unique combination of such zeta functions 
that leads to a variational formula like (|23"1), i.e. local up to cohomological terms. 

The variational formula for analytic torsion we obtained is more neatly expressed at the 
level of the Ray-Singer metric, since then the global term disappears. The next result is 
analogous to the variational formula for the Ray-Singer metric on Riemannian manifolds 



(see U, Theorem 4.14] and [|Tl|, Theorem 1.18]) 



Corollary 3.7. Let \\ \\c denote the contact Ray-Singer metric on det H*(£,dH)- 
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(1) The following identity holds: 

n „ 

(24) (In || || )- = V(-l) fc / tr(aa t p. fc )0A(d0r. 

(2) Under conformal variations of the contact form (6 e = e 2eT 6, J £ = J), for a function 
T, we have 

n . 

(In || \\c) 9 =2Y J (-l) k (n+l-k) / Ytr(a t o. fe ) A (d#) n . 
fc=o 

Proof. Recalling from fll8|) and §g the definition of the Ray-Singer metric, we have that 

(25) In || ||^ = 21nT c + ln| || 2(£) , 

where the L 2 metric | U 2 (£) is induced on det H*(£, du) from the inner product on H*(£, dn) 
defined by 



{[u],[v])e,j= / VuA*Tv. 
Jm 

But for the orthogonal projection V onto harmonic forms H*(£, djj), one checks that V* takes 
H*(£,dH) to its orthogonal complement. Thus 



(26) (M,M)*=/ VuA(*')Vv = ([u],a[v}). 

If we use Hodge * duality in the definition of detH*(£ ,dn), then take an orthonormal basis 
of each 7i k (£, dn), k = 0, . . . , n, and finally use (|26|), it is easy to see that 



(In | | c )- = ^(-l) fc Tr(aP fc ). 
k=0 

This together with (|25| ) and Theorem |Oj completes the proof of (1). 

Assertion (2) follows immediately from (1), since for conformal variations it is straightfor- 
ward to check that on £ k 

a = *~V = 2(n + 1 - jfe)T Id . 

□ 



Note that setting T = 1 in Corollary 3/7 (2), i.e. performing a constant rescaling 6 t— > e 2e 9, 
yields 

(27) (|| \\cY = 2k(0). 

In particular, if the contact invariant re(0) / 0, then we could not hope for any invariance 
of Ray-Singer metric. Note that by definition 



<0)=£(-l) fc (n + l-fc)C(A*)(0) 

k=0 

n 

= Y^(-l) k (n + l-k) tr(a t o. fe )0A(d0) n 

JM 



k=0 
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by Proposition 3.2, where again a t o. k is the constant t° coefficient in the diagonal small-time 



asymptotic expansion of the heat kernel of A k . Therefore n(0) is an integral over M of local 
curvature data, namely 

(28) k(0)= / P n (R,A,T, V)dvol, 



for some universal polynomial in Tanaka-Tanno-Webster ([g2[ 43, ^]) curvature i?, torsion 
A, Tanno's tensor T = V J, and their covariant derivatives. 



We show in Corollary 3^ below that, in dimension 3 (n = 1), k(0) vanishes identically. 
Whether contact invariants such as k(0) vanish in all dimensions is an open problem. For 
further discussion in the contact case see §7], and [0, Remark 9.3] for a similar problem 
arising for the eta function of the contact complex. 

3.3. CR/contact invariants in dimension 3. In dimension 3, besides the vanishing of k(0) 
we mentioned, we can also obtain more explicit variational formulae, and get CR/contact- 
invariant corrections to the contact Ray-Singer metric. 

Corollary 3.8. On 3- dimensional contact manifolds: 

(1) It holds that k(0) = 0, and thus || \\c is independent of a constant rescaling 9 i— ► K9 . 

(2) There exist universal constants C\, C 2 such that under a conformal variation (9 £ = 
e 2eT 9, J £ = J) we have 

(In || \\c) m = [ T(C 1 A H R + C 2 lmA u 11 )9Ad9, 

JM 

where R, A, An and a comma subscript denote respectively the scalar curvature, tor- 
sion, sub-Laplacian and covariant differentiation with respect to the Tanaka-Webster 
connection flg|, Hj °f (P, J )- 

(3) Let C[ = -C1/8 and C' 2 = C 2 /A, with Ci,C 2 as above. Then 

|| || C R = exp(c(/ R 2 9Ad9 + C' 2 [ \A\ 2 9Ad9)\\ \\ c 

\ JM JM ) 

is a CR-invariant (i.e. independent of contact form) metric on detH*(£,dH)- 

(4) There exist universal constants C3, Ca such that both 

|| ||fr = exp(C 3 i/(M))|| ||cR and || \\ H = exp(C 4 7?(^*))|| || C r 

are contact-invariant metrics, with v(M) the v-invariant of Biquard-Herzlich @, and 
rj(D*) the CR-invariant correction to the pseudohermitian eta invariant r](D*); see 
@, Theorem 9.4]. 

Proof. We complexify H and work in a local frame {Zi,Zj} and coframe j^ 1 ^ 1 }, where 
6 1 (T) = = 9 1 {T) (T the Reeb field). Under a constant scaling of contact form 9 = K9, the 
relevant heat coefficient a t o. k , for k = 0, 1, scales as tv(a t o. k ) = K~ 2 ti(a t o. k ). This is easily 
verified by an argument similar to that for (6.48)]. Basic invariant theory (see e.g. [41, 0) 
then tells us that tr(ct t o. fe ) must be a universal linear combination of 

(29) R 2 , \A\ 2 , A H R, R >0 = 2ReA u 11 and Imi u U . 



18 



MICHEL RUMIN AND NEIL SESHADRI 



Now k(0) is the integral of a linear combination of these terms, which is moreover indepen- 
dent of the choice of 8. A familiar argument (see e.g. j?]]) shows that the integral of a linear 
combination of R 2 and \A\ 2 can never be contact-invariant. Thus k(0) is the integral of a 
divergence, and hence vanishes. This, together with (^j), proves assertion (1). 

Consider now the differential of In || \\c under a conformal change of 9. This may be seen 
as a real 1-form a on the space of contact forms. By Corollary |3^ (2) and ( |29| ) it can be 
written 

a e (T)= [ r(ciR 2 + c 2 \A\ 2 + c 3 A H R + c±R + c 5 ImA 11 11 )6Ad6, 

JM 

for some universal constants Cj. Here we identified the tangent space TgQ with functions T 
on M. By [j?], Lemma 9.5], the general vanishing of such a 1-form on constant T implies that 
c\ = C2 = 0, while the fact that a is a closed form gives C4 = 0; see |@] for details. This proves 
assertion (2). 



Also by (83)-(84) in §9], one has 
_d_ 



/ R 2 0Ad8 = 8 [ T{A H R)9Ad6. 

JM JM 



[ \A\ 2 9Ade = -i[ TilmA^^eAdO, 

JM JM 



and 

_d_ 

dT 

leading to assertion (3). 

Assertion (4) is proved similarly as for the case of the contact eta invariant in |?], §9]. The 
CR deformations (i.e. of J) of the CR-invariants r](D*) and ln|| ||cr are all given by 
multiples of 

/ (Q,J)9Ad9, 

JM 

where Q is Cartan's tensor; see |@, §9] for details. □ 
Remark 3.9. As may be seen from (|24|), in order to determinate the various universal constants 



in Corollary 3.8 and investigate whether || \\c has any contact-invariant properties, one needs 
to calculate the local coefficients of t° in the diagonal small-time asymptotic expansion of 
the heat kernels of the fourth-order Laplacians we consider here. Formulae for calculating 
these coefficients are built into the pseudodifferential construction of the heat kernel, however 
implementing these in practice seems difficult. 

4. Contact and Ray-Singer analytic torsions of CR Seifert manifolds 

We follow @ to review the definition of a CR Seifert manifold and to fix notation. Note 
that in dimension 3 a calibrated almost complex structure J for the contact structure H is 
automatically integrable; the pair (H, J) is often called a pseudoconvex CR structure. 

Definition 4.1. A CR Seifert manifold is a 3-dimensional compact manifold M endowed 
with a pseudoconvex CR structure (H, J) and a Seifert structure ip : S 1 — > Diff (M) that 
are compatible in the following sense: the circle action cp preserves the CR structure and is 
generated by a Reeb field T. 
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It is easily proved that existence of a Reeb field T satisfying ip^{d/dt) = T is equivalent to 
existence of a locally free action of S 1 whose (never vanishing) infinitesimal generator preserves 
(H, J) and is transverse everywhere to H . 

The quotient space E = M/S 1 is an orbifold surface with conical singularities. Each CR 
Seifert manifold is then the S 1 -bundle inside a line orbifold bundle L over the compact Rie- 
mannian orbifold E. Singularities of L are located above the singularities of E in such a way 
that the total space M of the bundle is a smooth manifold: if the local fundamental group at 

■ 2-7T 

p £ E is Z/aZ (a € W), a generator acts on a local chart around p as e l ~ and on the fibre 

■ 2nf3 

above p as e « , where a and f3 are relatively prime integers with 1 < a < p. 

Theorem 4.2. Let M be CR Seifert manifold and p : tti(M) — > U(N) a unitary representa- 
tion. Then: 

• The analytic torsion Tc(p) of the twisted contact complex and Ray-Singer analytic torsion 
Trs{p) satisfy 

T c {p) = {T RS {p)Y l . 

• The two Ray-Singer metrics on det H*(M, p), corresponding to the de Rham and contact 
complexes (see (Q) and (|l9|) ), coincide, i.e. 

lie = II H-RS , 



via the isomorphism det H*(£* dn) = det H* (f2*M, d) coming from Proposition 2.1 



The remainder of this section will be devoted to the proof of these results. 

We first need to compare the two spectra coming from the de Rham and contact complexes. 
This has been done in [0, §§7, 8] in the untwisted case, i.e. for a trivial representation. We will 
rely on and refer to the spectral analysis done there and point out the few differences coming 
from the use of the representation p here. 

4.1. Circle action and Fourier analysis. Let V be the flat complex vector bundle over 
M associated to p. It is the quotient of the trivial bundle M x over the universal cover 
M of M by the deck transformations ^y.(m,v) = {r{^)m, p(^f)v). In what follows the contact 
complex is twisted by p in order to take values in V. 

Let (ft be the circle action on M induced by the Reeb field T. It may be lifted on V, by 
parallel transport for the flat connection V p , but no longer as a circle action. We have instead 

where / = <^[o,27r](^) is the generic closed orbit of the action, as seen in tti(M). This / is 
central, as comes from the presentation of the fundamental group of Seifert manifolds. 

Proposition 4.3 (see e.g. [20, |2^, |3£|). Let M be the circle V -bundle L of rational degree 
d = b + Yli % over the orbifold surface E of integral genus g, with n conical points xi of type 
(oti, f3i). Then ~n\{M) admits the presentation 

tti(M) = (/, a h bj (1 < j < g), gi (l<i<n)\ 

/] = hJ) = \9i, f] = g?T Pi = f b \{[a 3 , b 3 ] 9l = 1). 
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We can split V into irreducible components, where we have 

p(f) = e 2mx 

for some x £ [0, 1[. We recover a circle action on each such component V x by setting 

Using the circle action ip t one can still perform a Fourier decomposition as in |[|, §7]. Namely, 
for / £ V x , one has 

1 /* 27r 

/ = 5>»/ where vr n / = — / e~ mt Mf)dt 
nez 71 ^° 

satisfies tpti^nf) = e 7r n / and T(ir n f) = i(n + x)ir n f. Hence the spectrum of becomes 
the shifted Z, — x on V x . For A = — n — x we note 

(30) y A = v£ = v n {;r = -« - x} = 7r n (v x ) . 

As the circle action preserves the metric and the whole pseudohermitian structure (H, 9, J), 
we can split both the Hodge-de Rham and the contact complex spectra into their Fourier 
components. This is useful for comparing the spectra. 

4.2. Comparing the Riemannian and sub-Riemannian spectra. Following Proposi- 
tions 7.2 and 7.4 in @, we consider the following spaces. 

Definition 4.4. • Let TLy be the space of vertical 2-forms a = 9 A /3, with values in V, such 
that both a and Ja are closed. 

• Let also Tty be the space of pluri-CR functions in V, i.e. 

H v = kev(A H - T 2 ) = ker U V U V 

with Oy = d v dy and Dy = d* v dy p . 

According to ||, (63)-(68)] the non-zero spectrum of the non-positive second-order Laplacian 
P = D * +5hcIh on 2-forms splits into 

spec* (P) = spec* (£>*) U spec*(A#) 

(31) =spec*(A H ) U spec*(-Atf) 

Jspec*(-JT | n 2 v ) \ spec*(-|r| I 'Hy), 

where Ah = cIhSh + &h<Ih is the horizontal second-order Laplacian on functions. 

The torsion function kc of the contact complex is defined using the fourth-order Laplacians 
A 2 H on functions, and Ai = D*D + (cIhSh) 2 on 1-forms, with Ai conjugated to P 2 by Hodge 
* duality. Hence ([H]) yields 

(32) spec*(Ai) = 2 x spec* (A^) (J spec* (-T 2 \ Hy) \spec*(-T 2 | H v ) . 
Finally by ( |20| ) the torsion function of the contact complex reads 

kc(s) = 2((A 2 h )( S )-C(A 1 )(s) 

= C{~T 2 | H°v)(s) - C(~T 2 | H v )(s) + k(M, p ), 
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where we have set 

k(M, p) = 2dim(ker Ajj) — dim(ker Ai) 

(34) 

v ' = 2 dim(jy° (M,p)) - dim(H\M,p)) , 

since the twisted contact complex is a resolution computing the cohomology of M with values 
in V. 

We proceed similarly for the Hodge-de Rham spectrum, and work again with the calibrated 
metric g = d0(-, J-) + 6 2 . Set 



(35) Q ± = ±~ + ] Jj+A 

where Ao is Hodge-de Rham Laplacian acting on functions. According to |7|, Corollary 7.6] 
the spectrum of d* on 2-forms splits as 

(36) spec* (d*) = spec* (Q + )|J spec* (-Q~ \ {Hy) L ) (J spec*(- JT \ H 2 V ) . 
By Definition 0, we have A = A H - T 2 = \T\ - T 2 on H% so that 

Q- = -1/2+ VV4 + A = \T\ on Hy . 

Then (|^) also reads 

spec*(d*) = spec*((5 + ) spec*(— Q~) \ spec*(— \T\ \ Hy) spec*(— JT \ Hy) , 
and since 5d on 1-forms is * conjugated to {d*) 2 on 2-forms, we get 

(37) spec*(<5d) = spec* {Q + f (J spec* (Q~) 2 [j spec* (-T 2 | 7^) \ spec*(-T 2 | Ky). 
Now following our convention on analytic torsion, inverse to the original definition of Ray- 



Singer |33|] (see Remark \2J^ and (|i7|)), the torsion function of de Rham complex reads in 
dimension 3 as 

3 

krs(s) = ^(-l) fc+1 A;C(A fc )( S ) 

k=0 

= 3C(A )( S )-C(Ai)( S ), 
with Ai = d5 + <5(f the Hodge-de Rham Laplacians on i-forms. Using 

C(A 1 )(s) = C(5d)(s)+C(Ao)(s) 

and ( |37j ) one finds that 

K RS (s) = 2C*(A )(s) - C(8d)(s) +3dim(ker A ) - dim(kerAi) 
= 2C(A )( S )-C(Q + )(2 S )-C(Q")(2 S ) 

- C(-T 2 | + C(~T 2 | + k(M, p) , 

since kerQ + = {0} and kerQ - = ker Ao by the definition ( |35| ) of Q^. Comparing to the 
contact-complex torsion (|34]) we have shown the following result. 
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Proposition 4.5. On a CR Seifert manifold, the Ray-Singer and contact complex torsion 
functions twisted by a unitary representation satisfy 

(38) k rs (s) - k c {s) = 2C(Ao)(a) - C(Q + )(2s) - C(Q~)(2s) , 



where A is the Hodge-de Rham Laplacian on functions and Q = ±1/2 + a/1/4 + A . 
Note at this stage that the right-hand side of (|38|) vanishes at s = 0, as needed by the 



vanishing of both torsion functions at s = 0; see |33|] and Corollary ^J. This also follows from 
C(Ao)(0) = 0, for the Hodge-de Rham Laplacian in odd dimension, and that ((Q + )(0) = 
-C(Q~)(0), by |, Lemma 8.5]. 

4.3. Proof of Theorem |4.2|. In view of Proposition ^3] we need to show that if we set 

Q(s) = C(Q + )(s) + C(Q~) 00 " 2C(A)(s/2), 

writing A instead of Ao (the Hodge-de Rham Laplacian on functions) for brevity, then Q'(0) = 
0. We have a hint that this is true by examining finite energy cut-offs: at any finite spectral 
level (A < A) it holds that 

C(Q + )'(0) +C'(Q")'(0) = -lndet(Q+) -lndet(Q-) 

= -lndet(Q+ x Q~) 
= -lndet A 
= C(A)'(0). 

Hence Q'(0) is insensitive to finite eigenvalues and behaves like a pseudodifferential invariant. 
It may indeed be seen as a multiplicative anomaly for the regularized determinant of the 
product of the two commuting operators . As thus Q'(0) is related to a Wodzicky residue- 
type invariant; see |24, §6.5]. 

In fact the spectral function Q makes sense on any compact Riemannian manifold. The 
following result holds in a far more general setting than CR Seifert manifolds. 

Lemma 4.6. On any odd- dimensional compact Riemannian manifold Q(0) = Q'(0) = 0. 

Proof. Consider the one-parameter deformation 



Qf = ±A + V A 2 + A , 

so that, with A = 1/2, Qy 2 coincides with our original Q^. Note that the product formula 

Q+ x Q- = A 

we already mentioned is preserved during the deformation. By ellipticity of and a Mellin 
transform 

C(Qf)(s) = f^l t-'CTrCe-*^) -dimkerQ±)^ 
is holomorphic for large s. Define a function, holomorphic for large s, 



r+oo 

F±(A, S )= / ^Tr^e^A)^, 
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where here and in the sequel Tr*(P) = Tr(P) — P(const. function = 1). In particular 

i r+oo , 

C(^ /2 )( S ) = C(Ql /2 )(s) = f^y / o t s - lr iv(e-* Q v 2 )^ 

= i + r( s )- 1 F+(i/2, s ), 

and 

i r+oo _ 

C(Q- /2 )(s) = 1 + C*(Q V2 )( S ) = 1 + f^y J Q ^ Tr( e -* Q v 2 _ i)dt 

= i + r( s )- 1 F-(i/2, s ). 

Thus 

Q(s) = l + T(s)- 1 F + (l/2,s) + l + T- 1 (s)F-(l/2,s)-2-2T(s)- 1 F(0,s) 
~ s(F + (l/2, s) + F-(l/2, a) - 2F(0, a)) 
when s — > 0, with F(0, s) = F + (0, s) = F~(0, s). Therefore we need to show that 

F+(l/2,0) + F- (1/2,0) -2^(0,0) =0, 
for which it clearly suffices to show that 
(39) d x F + (X,0) + d x F-(X,0) = 0. 

Now one has, for the smooth family of commuting elliptic first-order operators Q\, 

■(e-*tf) =-t(±l + 



4(e^) = -t(±l + ; A )e-^ 



J " d ;( e -^), 



VA 2 + Adt 

so that 

/•+oo j /p-tQt — p-tQx \ 

or after integrating by parts, 

s / ^-iTr*f- 6 )rff 



= s 



/•-hoc /e~ ty A - e _ty A \ 

/ t'-i-iW* — ), 

/•+oo -tv/A^+A 

or after again integrating by parts, 

(40) d x F + {X,s) + d x F-(X,s) = 2s / g{X,t, s)Tr* (e' ty/ ^+ E )dt 

Jo 



with 

lu . 



g(X, t,s) = u s 1 smh(uX)di 
Jo 
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We therefore need to study the residue at s = of the integral expression in ( |40| ) . First g 
is easily expanded as 

q(X,t,s)= / u s_1 > — 77-du 

Jo ^ (2p+l)! 

_ A 2p+l t 2 P +l+ S 

(41) -^ o (2p + l)l(2p + s + l)- 

Consider next the Poisson kernel Tr*(e~*^ A2+A ) in (flO|); the beginning of its asymptotic 
expansion as t \ is related to that of the heat kernel Tr(e~^ A +A )) as follows. Recall from 
e.g. Lemma 1.7.4] that the trace of the heat kernel of a second-order elliptic Laplacian 
such as P = A 2 + A develops when t \ as 

fc>0 

where m is the manifold dimension and the Ck are integrals of curvature terms. 
Proposition 4.7. One has for P and Ck as above, for odd and even dimension m, 

[W 2 ] om-2fc rn _<)k + -\ 

(42) Tr*(e-^ 1/2 ) = £ L_ r( m 2^ + 1 ^^ _ 1 + /(t) 

where f[t)—*0 when t \ 0. 



,/vr 

fc=o 



Proof. This is a particular case of H, Theorem 3.1]. Indeed Bar and Moroianu gave there the 
full development of such kernels on the diagonal. Higher order terms in the development of the 
Poisson kernel are more involved in odd dimension since they contain log and even non-local 
coefficients. 

Here is an alternative proof of the partial development we need. The classical Laplace 
transform £(t _1 ^ 2 e _1// *) = \/irp~^/ 2 e~ 2p 1 leads to the subordination formula 

/•+00 



o 

between Poisson and heat kernels. Therefore summing at the trace level, 

f+OO 

iu 



r+oo 

Tr(e-* pl/2 ) = 7T-V2 / e -u u -i/2 Tl{e -t*P/4u )d 
Jo 

= 7T- 1 / 2 / e- u u- 1 ' 2 ( V c k t 2k - m (Au) m / 2 ~ k + B(t 2 /4u)) du 
Jo k=0 

l — n V JO 



k=0 



with B(t /Au) bounded and B{v) — > when v \ 0, This gives ( p2j ) by dominated convergence 
and the remark that 

Tr*(e-* pl/2 ) = Tr(e-* pl/2 ) - e~' A = Tr(e^ pl/2 ) - 1 + (1). 
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□ 

We can now complete the proof of Lemma 4J3. We split ( [Iq|) into 



d x F + (\,s)+d x F-(\,s) = 2s(j^ + J )g(\,t,s)Tr*(e- tpl/2 )dt, 



By (41) the second integral here is meromorphic with simple poles at s = —2n — 1 for n G N; 



in particular it is regular at s = 0. Set 



2 m-2k , m -2k + l\ 
C k = o ) C k 



7T 

so that by @ 

/ g(X,u,s)Tr*(e- upl/ )du = / V Y -. A rc^ 2 *""** 

Jo Jo ^ ^(2p+l 2p + l + s fc 

1/0 1/0 0<fc<[m/2]p>0 V F ; 

+ / g(X,u,s)(f(u) - l)du 







E E 



u , |(((J> (2 P +l)!(2p+l + s )(2p + 2 + 2t -,„ + ») 



+ holomorphic terms for Re(s) > — 1 . 



This expression has no pole at s = if m is odd, giving ( p9| ) and hence Lemma 14^. □ 



We now prove Theorem |42j. First Proposition [L5| , Lemma and ( |17| ) show that 

T c (p) = exp( K / c (0)/2) = exp(^ 5 (0)/2) = (T^p))- 1 . 

The equality of Ray-Singer metrics now comes from ( |l~8| ) and ( |19[ ) , using the equality of L 2 
metrics on H*(M,p) when H*(M,p) is represented by harmonic forms in the de Rham and 
contact complexes. Indeed these latter two notions coincide on CR Seifert manifolds because 



of vanishing Tanaka- Webster torsion; see [37, Proposition 12] 



5. The torsion function of CR Seifert manifolds and its dynamical aspects 

As an illustration of our viewpoint on analytic torsion, we first show how to compute it on 
any CR Seifert 3-manifold M equipped with a unitary representation p : 7Ti(M) — > U(N). 

As we will be only concerned with the contact torsion function kq in the sequel, we will 
denote it by tz for brevity. 

Surprisingly, on CR Seifert manifolds, we will see that the whole contact torsion function 
tz, not only k'(0), is expressible using topological data and combinations of Riemann-Hurwitz 
zeta functions, parametrised by dynamical properties of p with respect to the circle action on 



M. This leads in particular to a Lefschetz-type formula for the torsion; see Theorem 5.7. This 



extends a result obtained by Fried JrTj in the acyclic case using topological methods. 

In fact it turns out that our spectral torsion function k(s) may also be seen as a purely 
dynamical zeta function, constructed from holonomies along all closed orbits of the Reeb field 



T and its length spectrum. This will be shown in §§5.S and 5.4. 
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5.1. Torsion function k and the Riemann-Roch-Kawasaki formula. We first relate n 
to basic holomorphic data. 

Let V be the conjugate complex vector space to V, i.e. the same underlying real space with 
the opposite complex structure, and set 

W = V@V . 

Consider the dw complex 

(43) d w : Q°W = T(M, W) — ► Q^W = n^H ®W = T(M, A ' 1 ^* ® W) 
given on the CR (holomorphic) bundle W and let 

(44) 7i w = ker dw and 7i w = ker d* w 
denote its cohomology. Using — = iT, see e.g. 0, (57)], one gets 

spec(— \T\ | TCy) = spec(iT | H w ) , 

while, using Hodge * duality, 

spec(— JT | T~iy) = spec(iT | 7i w ). 
Then the spectral decomposition ( |3l] ) reads as follows. 

Proposition 5.1. The spectrum of P = D * +5hoIh splits as 
spec* (P) = spec* (A H ) U spec*(-A H ) 

[J spec* (iT | H w ) \ spec* 

where W = V © V and H w is the cohomology of dw as in (|43|) -([44]) . 

Remark 5.2. Compared to the trivial representation case treated in [j?], (68)], the only change 
here is the tensorisation by W. 

Now using the dw complex, equation fl32^ ) yields 

spec*(A 1 ) = 2 x spec* (A 2 H ) \J spec* (-T 2 | H w ) \ spec*(-T 2 | H w ) , 

and (|3^) becomes 

«( S ) = 2C(A 2 7 )( S )-C(A 1 )( S ) 
1 ° j = C(~T 2 | H^)W " C*(-T 2 | H w )(s) + k(M, p ) . 

This Lefschetz-type formula for k can be seen more topologically. Indeed, Fourier decompose 
each V x and let 

W= W x with W\ = V\®V\. 

Aespec(iT) 
Then using the holomorphic genus 

Xa(Wx) = dimH Wx ~ dimH Wx , 
the torsion function also reads as the Dirichlet series 

(46) k(s)= ^^ + *(M,p), 

ASspec* (iT) 
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where, from (|30|), spec(iT) splits into copies of — (Z + x) on each sub-representation V x of V, 
on which p(f) = e 2lnx . 

Remark 5.3. For comparison, the eta function of P = D * +5h<1h twisted by p may be 
expressed using Proposition 5^1 in a similar manner. One gets 

(47) „ P (P) W = Y. BignW Sg^jag) 

Aespec*(iT) ' ' 

which is strikingly the 'odd version' of the formula (46) for the torsion function k. Note that 
by 0, Theorem 8.8], 7] p (P)(0) identifies with r] (M,p), the diabatic limit of the Riemannian 
eta invariant with value in p 7 i.e., the constant term in the development of rj(M,g £ ,p) for the 
diabatic metrics (which we also consider in the present paper) g £ = e~ 1 d9 + e~ 2 2 . 

In order to express the series ( f46| ) using the Riemann-Roch-Kawasaki formula, we need 
first to see the bundles V\, a priori defined over M, as F-bundles over the orbifold S, and 
compute their degrees and orbifold exponents. 

Recall that we consider in ( |30| ) the circle action given by 

^ t = e- Ux vt = e int on V x = V? , 

where (fit acts by parallel transport along T in the flat V. Now e lt is a local coordinate on the 
fibre of the circle F-bundle M = S(L) — > 3D, hence e m * is a coordinate for S(L n ). Thus we 
have 

(48) V x = S(L- n ) ® F X . 

Sections of V^f are invariant by the circle action ip t and satisfy tpt = e ttx . Then one can see 
Vq as defined over S, but equipped with the non-flat connection 

for which sections s of V(f are parallel along T. 

Therefore, as seen from S, (V^jV 2 ) has curvature Q = —ixdO and rational degree (see 

ISi) 



<W) = ±^T±v*(-ixd6) 
= - dim{V x )xd(L) , 



because Ql = idO. Recall that by Proposition |4.3| , d(L) = b + Yli f 1 - Hence by (48) we get 
(49) d(V x ) = dim(V x ) (-nd(L) - xd(L)) = dim{V x ) Xd{L). 

By conjugation we have also 

d(\\) = -dim(V x )Xd(L) = -d(V x ) , 

so that finally 

d(W x ) = d(V^x) + d(V p , x ) = 0, 
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as expected for this smooth part due to the real structure on W\ = V\ @ V\; see |26|, §14]. We 
did the above computation for completeness, as the degree of V\ is needed to study the eta 
function given in Remark |5.3| . 



We determine now the local action on V\ of Tj = Z/c^Z around an orbifold point G S. 
Recall that locally over Xi, M = S(L) is the quotient of C x S 1 by the action of Tj generated 
by 

g.(z 1 ,z 2 ) = (e 2 ^z 1 ,e 2 ^^z 2 ). 
Since z 2 = e lt here, we get the local action on VJf 

g.(z x ,z v ) = (e^/^zt^Ttfo/a^zv). 
Let fi = (f([0, 2n/a%]) be the exceptional closed primitive orbit over Xj. We have 

(50) =e~ i ' 1 '«^ x l ai p{$if\ 

As = / in 7Ti(M), p(fi) ai = p(f) and the spectra of p(fi) satisfy 

( spec p(fi) = {e 2i ™^ | 1 < j < dirndl 

(51) < x + ki ■ 

I with Xij = — G [0, 1) and kij G Z . 

v on 

Hence ( |50| ) means that spec(g) = { e 2Mrfc ijAM} on the fibres of Vq around Xi, and, by ten- 
sorisation with L~ n in (|48|), the isotropy exponents of the action on V\ are all the couples 

(ai, (—n + kij)fa mod (Xi). 

Now the Riemann-Roch-Kawasaki formula (see [^0], [2^, [?]]) states that 
Xs(W x ) = dim(W A )(l - 9 ) + deg(W A ) - ^{^|} 

(52) = dim(^) X (S) - W + + / ("-^ 

where {a} = a — [a] G [0, 1) is the fractional part of a and = 2 — 2g is the Euler 

characteristic of the smooth surface £ associated to the orbifold S; see e.g. |29|. Observe that 
(|52]) does give integers since {a} + {—a} is when a G Z and 1 otherwise. Recall also that, to 
ensure smoothness of the ^-bundle M = S(L), the numbers a, and fa are assumed relatively 
prime, thus giving a free action of Z/a^Z at orbifold points. Hence the fractional part in ( |d2"| ) 
simplifies using 



(53) 8(n,i,j) 



1 if n — kij G (XiL 
otherwise . 



Then we have 

(54) Xa(Wx) = dim(F*)x(E*) + £<5(n,t,i) 
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where 

X(£*) = 2 - 2g - \I\ 

is the Euler characteristic of the punctured surface X* = £\U/{xj} at the \I\ orbifold points. 

For a s]0, 1] let ((s,a) = — be Hurwitz zeta function. We can now express k as 

' (n + a) s 

a combination of such functions. This is the first step towards the identification of the torsion 
function dynamical zeta function given in § |5,3| . 

Theorem 5.4. Split V into irreducible V x ; then the torsion function spectrally decomposes as 



E* 

such that : 

• On V x with x €]0, 1[, i.e. p(f) = e 2inx + Id, we have 

k x (s) = dim(V*) X (E*) {((2s, x) + C(2s, 1 - x)) 
(55) + E (C(2s, x^) + C(2s, 1 - xy)) . 

• On V° = ker(Id - let = ker(Id - p(fi)); then we have 

«o(s) = «(M,p)(2C(2 S ) + 1) + 2C(2 S )^dim(F ' i )(a- 2s - l) 

(56) 



+ E -^{C( 2s ^i,j) + C^s,l-Xij)) . 



a, 

This relates the torsion function to dynamical properties of the circle action here. Indeed 
apart from the cohomological term n(M,p), the expression is clearly built on the holonomy 
properties of p along the various closed primitive orbits of the flow: the generic orbit / of the 
action over £*, and associated holonomy p(f) = e 2l7TX on V x , and the exceptional orbits fa of 
holonomy p(fc) = {e 2wxi 'i}. 

Proof. • We compute first the contribution of V x for x ^ 0, i.e. when p(f) = e 2l7TX / Id. Here 
iT = A = — n — x / always, and by @) and 



= dim(^) X (S*) £ — L_ + E E 



1 



,n + x| 2s ^— ' ^— ' Ifca,- + fc< i + x\ 2s 
nez 1 1 i,j fcez 1 1 lJ 1 



In + x| 2s ^— ' I — n + x| 2s 

n>0 1 1 n>0 1 1 



+ ^ a 2s S lit + xijl 25 ' 



by (51). This leads to (fa). 
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• We compute now Kq, including the cohomological term k(M, p) from ( |34| ) and (46), since 
harmonic forms only appear in ker(iT) C V P: o; see e.g. |37|, Proposition 12]. We have 

(57) *„(,) = dim(y°) x (s*) E ri + E E 6J j$f- + *( M < p) ■ 



We recall from (|3J) that 

k(M, p) = 2 dim H° (M, p) - dim H 1 (M, p) 

can be computed using contact-harmonic forms on M. By |37], Proposition 12] contact- 
harmonic forms are both holomorphic and T-invariant since the Reeb flow preserves J here, 
i.e. Tanaka-Webster torsion vanishes. Then one gets 

k(M, P ) =2xq(W ) 

= dim(y°) X (S*) + ^5(0,i,j) by 

hi 

(58) k(M, p) = dim(y°)x(S*) + ^ dim(y°' i ) , 



since, by (pl|) and (p3|), 5(n,i,j) = 1 if and only if Xij = 0. Then (|57|) reads 



o,-. 



Ko (s) = K(M,p)(2C(2 S ) + l) + ^ ^ ^^l_2C(2 S )^dim(y 

i,j neZ* i 

We observe now that if Xij £ ]0, 1[, 

S(n,i,j) v-^ 1 



lnl 2s / / 



^ I (C(2 S ,x ij ) + C(2s,l-x i , i )) by (ID. 



On the other hand for Xij = 0, 



\ - 6(n,i,j) _ \- 1 _ J2_ v 
^ \n\ 2s ^ \kai\ 2s a 2s ^ >' 



as needed in (56). □ 

The expression of k x given in Theorem vanishes at s = 0, as it ought to by Corollary 
3.8 . Here this follows from the classical result C(0>°) = 1/2 — a; see |^6|, §13] for instance. The 
following observation will be useful in the sequel. 

Corollary 5.5. The torsion function k has a unique simple pole at s = 1/2 with residue 

Res 1/2 (K) = x(£)dimV, 
where x(E) = 2 — 2g + $^(— — 1) denotes the rational Euler class of the orbifold E. 
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Proof. One knows (see J46], §13]) that the Hurwitz zeta function ((2s, a) has a unique simple 
pole at s = 1/2 with residue 1/2. Then ( |55|) on V x yields 

Res 1/2 (K X .) = dimV x X (Z*) + dimF E ^— = dim V X X {X) . 

On the other hand by ( |56| ) in F , 

Res 1/2 ( K0 ) = n(M,p) + Vdim^'^C- - 1) + V dim((y°' i ) ± )- 

= (k(m, p)-J2 dim y0,i ) + dim v ° Yl — 

i i 
= dimy°( X (S*) + Y, —) = dimy° X (S) , 



I 



by <m. □ 



Remark 5.6. We lastly observe that a similar treatment applies to handle the twisted eta 
function in Remark |5.3|. Indeed by the Riemann-Roch-Kawasaki formula and (E9|) one has 



Xa(Vx) - Xa(Vx) = 2dim(V x )Xd(L) + £ 



n + kij)(3i \ f (n — hj)f3; 



on on 



and by ( fl7| ) the contribution of V x to eta is 

ifc(P)( a )=2dm(V*)d(L) £ ^ 

Aespec* (iT) ' ' 

+ E E ( 2 |(=-±*u)A|_ 1 + 4(Bi4ij) )^ 

ij Aespec* iT I- 1 J II 

The generic smooth contribution may be written as 



2dim(y :r )d(L) x 

taking value 



2((s - 1) if x = 

C(s - l,x) +C(s - 1,1 - x) if 



-dim(V*)d(L)(i + x(l - x)) 
6 

at s = 0; see |4|, §13]. Following Nicolaescu's work pEj] , the remaining 'periodic' eta term can 
be handled using Dedekind-Rademacher sums and Hurwitz functions; see Proposition 1.10 
and Lemma 1.11 in |29] for details. 
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5.2. A Lefschetz-type formula for the Ray— Singer metric. Using Theorem 4.2 we can 
now compute the Ray-Singer analytic torsion Trs = exp(— «'(0)/2) = (Tc)" 1 , which gives 
the associated Ray-Singer metric on det H*(M, p), 

(59) || Has = (Trs)' 1 \ \l 2 (q*m) ■ 



In the acyclic case, i.e. H*(M,p) = 0, Fried [|T7j] has shown that the Reidemeister-Franz 
torsion, and thus the analytic torsion by the works [15, 28] of Cheeger and Muller, may be 



nicely expressed 'a la Lefschetz' using determinants associated to the generic and exceptional 
holonomies along the primitive orbits of the circle action. For a general unitary representation 
p : 7Ti(M) -> U(N) we obtain: 

Theorem 5.7. Let p(f) T and p(fi) T denote the restriction of these holonomies to respectively 
(V ) 1 - and (V ' 1 ) 1 - with 

V° = ker(Id - p(/)) and V ^ = ker(Id - p(fi)) . 

Then Ray-Singer analytic torsion Trs(M, p) = exp(— k (0)/2) is given by 

(60) T RS (M,p) = ^^^-pm^U^^^^- 

Proof. By Lerch's formula d s C(s,x) s= Q = lnT(x) — |ln(27r), see |^6|, §13], we have 

d s ({0, x) + 9 S C(0, l-x) = ln(r(x)r(l - x)/2?r) 

= — ln(2 sin(7T3;)) by the Euler reflection formula, 
= -ln\l-e 2inx \. 

Hence by (||) on V x , 

-4(0)/2 = dim(y x ) X (S*)ln|l _ e 2 ^| +^ln|l -e 

hi 



2inxij I 



which gives the determinant contribution of V x to ( |60|) . 
By fl56| ) on F° and Lerch's formula again, one finds 

-k' (0)/2 = -2C'(0)/c(M,p) - ^dim(T/°'')lnai + ln|l-e 2 

i i,j\x iy j^0 

similarly as above. This gives the needed contribution of V° to (|SC|), □ 

As required, formula ( |60| ) coincides with that of Fried |l7[ p. 198] for acyclic representations. 
The only new factor in our case is the cohomological term (2iz) K ( M > p K That the full expression 
for the torsion is 'quantised' here is due to the rigidity of volume in this CR Seifert case. 
Namely, the size of 6 is fixed such that the circle action is generated by the Reeb field T in 
constant time 2tt, hence the volume forms dvol = 9 A d9 on TM and d9 on H are also fixed. 
Then when the complex structure J changes on the base E, together with the calibrated 
metric g = d6(-,J-) + 2 , harmonic forms representing H*(M,p) change, but their length in 
detW*(M, p) does not. Hence the 1? metric used in ( |59[ ) is independent of J, as is the analytic 
torsion given in (|60|). 
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5.3. The contact torsion function as a dynamical zeta function. In |l6], [l?]], Fried 
proposed to express the torsion using the following basic dynamical objects. 

For each free homotopical class C of periodic orbit of the Reeb field T, let 1(C) denote its 
length and ind(C) its Fuller index; see ||l9| or [JR], §4] for an account of these notions. 



Proposition 5.8 (||16, Lemma 5.3], |T7|])' The free homotopy classes of closed orbits ofT are 
the following : 

(1) f n with n 6 N* ; of length 2nn and Fuller index x(E)/n ; where 

X (E) = 2 - 2g - - l/cti) = X (S*) + £ l/a t 

i i 

is the rational Euler class of the quotient orbifold S = M/(T); 

(2) the isolated /" for n £ a^N, of length 27rn/aj and Fuller index 1/n. 

Fried observed that for acyclic unitary representations one has 
(61) Tfis(M, p) = |exp(Zp(0))| , 

where Zp(0) stands for the analytic continuation at s = of the dynamical function 

Z F ( S ) = -J2™d(C)Tt(p(C))e- se ^ . 

c 



This can be checked directly from the calculation of torsion in Theorem |5/7| and Proposition 
5.8 , as in [[l?], §1]. Such a link between analytic torsion and flow dynamics is not coincidental; 
it has already been observed in many other geometric situations, see e.g. [16|. In particular it 
holds for the geodesic flow on hyperbolic manifolds, as proved by Fried in [18], or more generally 
on locally symmetric spaces of non-positive sectional curvature, as proved by Moscovici and 
Stanton in |27|. These results are rooted in Selberg's trace formula, expressing heat kernel 
traces as a sum of traces along closed geodesies. 

In our setting, in view of Theorem 5.4 , it is quite natural to try to express the contact 
torsion function k(s) itself using the same dynamical data as above. Indeed the whole spectral 
function k may be nicely interpreted 'a la Selberg' as a purely dynamical zeta function of the 
Reeb flow. 

Theorem 5.9. Let 

/00 = r(«)co8( T ). 

For C closed let R Tr(p(C)) denote the real part of the trace of p(C) on V p . Then 

(62) f(s){ K (s/2)-K(M,p)) = ^ind(C7) fi Tr(p(C7))£(C) s , 

c 

where the sum is taken over all free homotopical classes of closed orbits of the Reeb flow. 

We first make some comments about this identity. First, to remove the real part in (|62}), one 
could also sum over all orbits C and C , the latter corresponding to the opposite flow — T. 
Indeed the torsion function is not sensitive to a change of 6 i— > —6, together with J i— ► — J, 
since k is defined using real operators Aq and Aj. 
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Also let 

(63) Z p {s) = md(C) R Tr(p(C))£(C) 

c 



be the dynamical zeta side of fl62|). Prom Proposition 5^ this converges for Re(s) < 0. In 
contrast the spectral side 

K*(a/2) = k{s/2) - k{M, P ) = 2Tt*(A~ s/2 ) - Tr*(A~ s/2 ) 

is a converging series for Re(s) > 2. Hence, when seen as series, the spectral and dynamical 
sides of (32) never converge for the same s, and the identity only holds through meromorphic 
continuation. When s — > we get 

(64) jm( Zp ( a ) + <^)) =k'(0)/2 = -HTrs(M, P )), 

and (In of) the analytic torsion may be seen as a topological regularisation of the formal 
dynamical series 

(65) -" Z p (0)" = -J2 md(C7) R Tr(p(C7)) . 

c 

Comparing with (|6l| ) and Fried's dynamical function yields 

"Z p (0)" = Z p (0) = -Re(Z F (0)), 

in the acyclic case, and the dynamical functions Zp and Z p both provide analytic continuation 
of the same dynamical series in (|65[). This series has been interpreted in [[0], [lTj as being the 
total Fuller measure of periodic orbits, and has a formal invariance by deformation of the flow, 
as long as orbit periods stay bounded. 

We note also that the trace formula ( |62|) can be written in a more symmetric manner 

(66) r(sK( S ) = ^T(~ - s)Z p (2s) , 

V 71 " ^ 

as follows from the classical identities (see p6]| ) 

T(s)T(s + h = 2 1 - 2s ^n^) and T(s + hr(-s + \) = . 

2 2 2 cos(7ts) 

This formulation will be useful in SI5.4. 



As a last comment, we observe that the trace formula (|62|) is homogeneous in the constant 
rescaling 9 \— * K6. Indeed, the metric here is g = d6(-, J-) + 6 2 , hence 1(C) = j c 6 changes to 
K£(C), while the fourth-order contact Laplacians Aj are homogeneous and rescale to K~ 2 Ai. 
Thus C*(A 1 ^ 2 )(s) rescales to K S C* (A l / 2 )(s) as needed. Such a property does not hold for the 
Hodge-de Rham Laplacians. 

Proof of Theorem [5j]. We will start from the expression of n(s) by Hurwitz zeta functions as 
given in Theorem 5A. First Hurwitz's formula (see §13]) states that for Re(s) < 

2r(l — s) r /7rs N v^ ,2irxn^ is ^ ,2nxn^ 

C(s,x) = —. — n sin — )> cos — = + cos( — ) > sin — = 

sv ' 1 (2vr) 1 - s L y 2 J ^ l n 1 - s/ K 2'^ V n i- s ) 

n=l n=l 
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so that using /(s)/(l — s) = tt/2 gives 



(67) f(s)(((s,x) + ((s,l-x)) = J2M e 



2inxn\ (27Tn)' 



n 

n=l 



Note also the corresponding limit expression when x — > + with Re(s) < 0: 



(68) 2/( S )C(s) = E 



n 

n=l 



We study the contribution of k x on V x with (0,1). By (p7[), formula ( pq ) yields 
fM r W ^^V- RlV W n )) .o ^ i ^y^ R TV(p(/D) ,2^n^ 



n>l i n>l * 



hence by Proposition p\q , 

/(*K(s/2) = ^(ind(D - £ 1) R Tr(p(r)K(r ) 4 



n>l i 



+ E E md(/r) R Tv(p(/f)K(/f) 



i fc>l 

= ^ind(C) R Tr(p(c7))/(C) s , 

C 

as needed in j6^). 

• We now study k on V° = ker(Id - p(/)). Recall that specp(/j) = {e 2 "™^} and F '* 
ker(Id-p(/i)). By formula © reads 

(69) f(s)(K (s/2) - k(M,p)) =k(M, P ) — + 2^2^ ~ (— J 

n>l « n>l 



+EE dim (^)(«r-i) (2vrn)s 



n 

i n>l 



By (|58|) , the first series reads 

k(m, P ) e ^ = [dim(v°)(x(E) - e h + E dim ^)] E (27rn 



n * — ' «i L — ' * — ' n 

n>l i i n>l 



E md(D K Tr(p(f n )W n Y "EE dim(F )^^ 

n>l i n>l 



n 

i n>l 
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Since p(fi) = Id on V ' 1 , the second series in ( |69| ) splits into 

RTV^o.^x^/f)) 2vn s _ s-^s-^ R Tr(p(f?)) 2-Kn y dim V * 2im , 8 

i n>l i n>\ i n>\ 

= E E md(/r) R TY(p(/r))w) s 

j fc>l i n>l 1 

since ^(/J 1 ) = /o(/ fc ) = Id on V° for n = fcaj. Therefore after cancellations ( |69| ) yields 

/( S )(« (s/2) - k(M, P )) = ^ind(C) R Tr(p(C7)K(C7) s , 

c 

as needed. □ 



5.4. The torsion heat trace as a dynamical theta function. In Theorem |5.9| , spectral 
and dynamical aspects of analytic torsion are compared through zeta functions. One can also 
work at the level of heat kernels. Consider the heat operators of the fourth-order Laplacians 
Ao and Ai of the contact complex, and set 

(70) Tr K (e- tA ) = 2Tr(e-* Ao ) - Tr(e~* Al ) 

= Tr* K (e- tA )+ K (M,p). 

Recall that for Re(s) > 1, 

i r+oo 

(71) k*(s) = 2C(A )(s)-C(A 1 )(s) = — -r / t^TrKe-^dt. 

Then the following trace formula holds in our CR Seifert setting. 
Theorem 5.10. One has 

(72) Tr K (e-* A ) = dim + _L V t(C) ind(C) R Tr(p(C))e- e <& 2 /4* , 

y/t y/nt q 

where C runs over free homotopical classes of closed orbits of the Reeb flow and x(£) is the 
rational Euler class of the quotient orbifold. 

Hence the torsion heat trace may be compared with the purely dynamical theta function 

(73) 0(t) = --L V £(C) ind(C) R Tr(p(C7))e-^ c ) 2 / 4 ' . 

We will first need the following fact on the asymptotic heat development. 
Proposition 5.11. As t \ 0, it holds that 

Tr K (e-* A ) = dimy^S^ + O(Vt) . 

y/t 
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Proof. By Corollary |5.5| , the torsion function k has a single simple pole at s = 1/2 with residue 
x(S)dimF. On the other hand we know that, for a fourth-order hypoelliptic Laplacian in 
dimension 3, as t \ 0, 

Tr K (e-* A ) = | + ^+ C0 + O(^). 

Mellin's transform ([Fl]) split into f + f^°° providing 

ci = Res s= i(«(s)) = 0, c 1/2 = r(l/2)Resi/ 2 (K) = \/7rx( s ) dimV , 

and 

c = k(M, P ) +Res (r(s)K*(s)) = 0, 
since «*(0) = «(0) - k(M, p) = -k(M, p). □ 

We can now prove Theorem 5.1C| . 
Proof. One takes Mellin transforms M. of both sides in (|72|). For Re(s) < one finds 

(74) ^ W(s) = ^!r(^- S )Z p (2 S ). 

On the other hand (|7l| ) and Proposition [5.11 yield that, for Re(s) > 1/2, 

In order to compare these identities we need first to extend them to a common domain. Indeed 

set 

Tr (e-' A ) = Tv* K (e- tA ) + («(M,p) - dhnV^^)x m (t); 



by Proposition |5.1l| this has a Mellin transform for Re(s) > —1/2 and 
(75) M(Tr (e~ tA ))(s) =T(s)k*(s)+ k(M ' p) +dimF' / ~ 



s- 1/2 
Set also 

Mt) = m - (<M,p) - dimy ^ s) ) X[1 , +oo[ (t) , 



so that (74) yields, for Re(s) < 0, 



M{ms) _ ^ r( i - s)Zp (2 S ) + ^ + dim V--^< E) 



/tt 2 I I s s- 1/2 

M(Tr (e-* A ))( S ), 



for —1/2 < Re(s) < 0, by ( pq ) and (|74|). Hence by injectivity of the Mellin transform, coming 
from Fourier injectivity on integrable functions here, one concludes that $o(£) = Tro(e~* A ), 
yielding the trace formula. 

The authors thank Patrick Gerard for an enlightening discussion on this proof. □ 
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Formula ( f72| ) is a typical Selberg-type trace formula, which holds in many other geometric 
situations, see (ll| p. 57] for instance. Note that it holds here even in variable curvature. 
This may look unusual as Selberg's technique is algebraic and relies on group actions, hence 
makes sense on uniformised (locally symmetric) manifolds. However we know, by (|46| ) or 



Theorem Hj], that on CR Seifert manifolds the torsion function k(s) is 'topological', meaning 
independent of the complex structure J since 9 is fixed by the circle action. Now except 
for some cases that fibre over the sphere S 2 with two singular points, all other CR Seifert 
manifolds can be uniformised, i.e. endowed with a constant curvature metric; see e.g. |H| or 



|20, Theorems 1.1, 1.2]. By Moser's lemma this can be done with a fixed volume. Thus, except 
for the special cases mentioned, one could have worked over a uniformised orbifold S, where 
Selberg's technique should also be applicable. 



The trace formula ( |72j ) has a striking consequence for the small time behaviour of the torsion 
heat trace Tr K (e~ tA ). Indeed, from its definition (|73|) and Proposition 5J, the dynamical theta 
function $ clearly decays very fast as 

= 0{e~ c/t ) when t\0, 



so that instead of Proposition 5.11 we get the full heat development 



Corollary 5.12. On CR Seifert manifolds, as t \ we have 

Tr K (e~ tA ) = dim + 0( e - c /') . 

V* 

Thus on such manifolds this torsion heat trace does show a 'fantastic cancellations' phe- 
nomenon as encountered for the heat of Dirac operators. 

Note that the only surviving term as t \ in this heat development is the integral of 
curvature data, as should be the case. Gauss-Bonnet reads here 



/ R0Ad9 = 2n f Rd6 = 4^ 2 X (S) , 



where R stands for the Tanaka- Webster scalar curvature, which coincides with the Riemannian 
scalar curvature of the base in this Seifert case; see e.g. 0. Then, by universality of the 
coefficients of the heat development, the following 'fantastic factorisation' holds on any 3- 
dimensional contact manifold. 

Corollary 5.13. On any contact 3-manifold, the full development of Tr K (e~ tA ) as t \ is 
of type 

(76) Tt K (e- tA ) ~ [ R6Ad9 + Yt n / 2 [ P n (R, A)dvol , 

where all invariant curvature polynomials P n (R,A) factorise through Tanaka- Webster torsion 
A = CtJ and its covariant derivatives. 



In the opposite direction, when t — > +oo, the trace formula (|72j) gives the asymptotic 
development of the dynamical theta function Indeed, after removing the zero eigenspace, 
heat decays exponentially, yielding the following property 
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Corollary 5.14. On CR Seifert manifolds, it holds as t —* +00 that 

(77) 0(t) = s(M, p) - dim V ^ X ^ + 0(e~ ct ) . 

This behaviour is not surprising; it can also be seen from the explicit formula ( |7q) and 
Proposition |5.8| , which relate the dynamical theta function to the classical Jacobi theta func- 
tion, whose decay at +00 is well known; see e.g. j2|, §21.51]. 
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